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ABSTRACT 

We develop an analytical model to follow the cosmological evolution of magnetic 
fields in disk galaxies. Our assumption is that fields are amplified from a small seed field 
via magnetohydrodynamical (MHD) turbulence. We further assume that this process 
is fast compared to other relevant timescales, and occurs principally in the cold disk 
gas. We follow the turbulent energy density using the Shabala & Alexander (2009) 
galaxy formation and evolution model. Three processes are important to the turbulent 
energy budget: infall of cool gas onto the disk and supernova feedback increase the 
turbulence; while star formation removes gas and hence turbulent energy from the 
cold gas. Finally, we assume that field energy is continuously transferred from the 
incoherent random field into an ordered field by differential galactic rotation. 

Model predictions are compared with observations of local late type galaxies by 
Fitt & Alexander (1993) and Shabala et al. (2008). The model reproduces observed 
magnetic field strengths and luminosities in low and intermediate-mass galaxies. These 
quantities are overpredicted in the most massive hosts, suggesting that inclusion of gas 
ejection by powerful AGNs is necessary in order to quench gas cooling and reconcile 
the predicted and observed magnetic field strengths. 
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1 INTRODUCTION 

Radio synchrotron emission of high energy electrons in the 
interstellar medium (ISM) indicates the presence of mag- 
netic fields in galaxies. Rotation measures (RM) of back- 
ground polarized sources indicate two varieties of field: a 
random field, which is not coherent on scales larger than 
the turbulence of the ISM; and a spiral ordered field which 
exhibits large-scale coherence (e.g. Stepanov et al. 2008). 
For a typical galaxy these fields have strengths of a few 
yitG. In a galaxy such as M51, the coherent magnetic field 
is observed to be associated with the optical spiral arms 
UPatrikeev et al. 2006(1 . Such fields are important in star for- 
mation and the physics of cosmic rays, and could also have 
an effect on galaxy evolution, yet, despite their importance, 
questions about their origin, evolution and structure remain 
largely unsolved. 

The Square Kilometre Array (SKA) will help us answer 
questions such as these. The SKA will observe polarized syn- 



chrotron emission (Gacnsler et al. 2004), and the "All-Sky 
SKA Rotation Measure Survey" will expand RM data sets 
by five orders of magnitude, mapping the magnetic fields 
of galaxies in unprecedented detail ( |Stepanov et al. 2008| 
Gaensler 2006). In particular, the SKA will provide data 
on the evolution of galactic magnetic fields to high redshift 
IjGaensler 2 006). making a theoretical model for this pro- 
cess very valuable. Here we present such a model within the 
context of hierarchical structure formation. 

Our model is based on the observed properties of galac- 
tic magnetic fields. Observations and simulations suggest 
that the random field is generated by turbulence in the 
ISM, which is modeled as a single-phase magnetohydro- 
dynamic (MHD) fluid, within which magnetic field lines 
are frozen (e.g. Cho et al. 2009). Simulations have shown 
that in such a turbulent MHD fluid, the random magnetic 
field energy and turbulent fluid kinetic energy are approx- 
imately equal after a few turnover times (|Cho et al. 2009). 



2 Shabala, Mead & Alexander 



so by writing an equation for the turbulent energy in the 
ISM the random magnetic field energy can also be evalu- 
ated. The large-scale ordered field is produced by the dif- 
ferential rotation of the galaxy, which winds the random 
field into a spiral - this is the basic operation of a dy- 
namo. A simple model for this process is postulated. The 
energy sources which contribute or remove turbulence from 
the ISM are supernovae, star formation and accretion of gas 
from the hot gas halo, and these parameters can be de- 
termined from well-established semi-analytic galaxy mod- 
els, in which galaxies form from gas condensing at the cen- 
tres of hierarchically merging haloes ([White fc Rees 1 978 
|Lacey fc Silk 1991| ICroton et al. 2006[) . The formation of 
magnetic fields in galaxies is intimately linked to galaxy for- 
mation, since the same physical processes are at work in 
both cases. 

Recently, Arshakian et al. (2009) have presented a qual- 
itatively similar model for the evolution of magnetic fields in 
late-type galaxies. These authors considered three processes. 
At high redshift, the initial seed field of ~ 10 -18 Gauss is 
amplified via the Biermann battery mechanism to generate 
what they refer to as the regular field. This field is rapidly 
(on eddy turnover timescale, ~ 10 7 years for the Milky Way) 
amplified by virialization turbulence. The mean-field galac- 
tic dynamo mechanism then amplifies this regular seed field, 
with a typical e-folding time of ~ 10 8 years, until it reaches 
equilibrium with turbulent dissipation. By contrast, given 
the short (compared to cosmological) timescales involved, 
in our model we simply assume instantaneous equality be- 
tween the magnetic field energy and turbulent kinetic energy 
to estimate what we refer to as the random field. This field 
is then ordered by the mean-scale dynamo on much longer 
timescales. 

While our treatment of large-scale dynamo generation 
is simplified compared to that of Arshakian et al., it encap- 
sulates the relevant physics and has the advantage of being 
simple enough to be included in a cosmological framework. 
In that sense this work is highly complementary to the Ar- 
shakian et al. (2009) study. We present here the first (to 
our knowledge) attempt to include magnetic fields in a self- 
consistent galaxy formation and evolution model. A number 
of galaxy properties are predicted, and we compare these 
with available data. This paper primarily focuses on the ra- 
dio properties of local late-type galaxies. 

The paper is organised as follows. In Section [5] we out- 
line the adopted galaxy formation model, and develop mod- 
els for both the random and ordered magnetic fields. Sec- 
tion investigates the salient features of our model and 
constrains various parameters. Predictions of our model are 
compared with available observational data in Sections [3] 
and[S] We summarise our findings in Section [5] 

Throughout the paper, we adopt a flat cosmology of 
Q,m = 0.3, Qa = 0.7, h = 0.7 and a s = 0.9, consis- 
tent with the 2dF GRS (|Colless et al. 2001|) and WMAP 
USeljak et al. 2005[ ) results. 

2 TRACING GALACTIC MAGNETIC FIELDS 
2.1 Galaxy evolution model 

We employ the analytical galaxy formation and evolution 
model of Shabala & Alexander (2009). This traces gas cool- 



ing, star formation, and various feedback processes in a cos- 
mological context. The model simultaneously reproduces the 
local galaxy properties, star formation history of the Uni- 
verse, the evolution of the stellar mass function to z ~ 1.5, 
and the early build-up of massive galaxies. While details can 
be found in Shabala & Alexander (2009), below we briefly 
outline those features of the model pertinent to the current 
paper. 

The model follows average properties of a galaxy of a 
given mass at redshift zero. Mass accretion history of the 
host halo is followed via the formalism introduced by van den 
Bosch (2002), providing analytical fits to detailed N-body 
simulations. Following the conventional picture, baryons are 
assumed to be accreted onto the halo together with the dark 
matter. Thus, in our model each halo with a given redshift- 
zero mass can only host one galaxy. Reionization feedback 
from first generations of structure limits the accreted bary- 
onic fraction. The accreted gas is shocked to the virial tem- 
perature and density, and cools radiatively. Once cold, the 
gas is added to the galaxy disk on a dynamical timescale, 
where star formation can take place. 

These cooling processes complete with various modes of 
feedback. Stars more massive than 8 Mq inject supernova 
shocks into the surrounding gas as they expire; this type of 
feedback is important in low and intermediate mass galax- 
ies. Gas cooling also fuels intermittent AGN activity, which 
becomes important in haloes more massive than 10 12 Mq. 

Such detailed treatment of the gas heating and cooling 
processes as well as star formation in an analytical frame- 
work makes this model ideal for following the cosmological 
evolution of galactic magnetic fields, as we outline below. 

2.2 Magnetic fields in galaxy disks 

Observations of spiral galaxies, for instance using faraday 
rotation measures, have indicated the presence of random 
and ordered magnetic fields, of magnitudes of a few fj,G (e.g. 
Stepanov et al. 2008). The random fields are coherent on 
small scales, of order parsecs. The ordered fields are coherent 
on much larger scales, and are observed to follow the spiral 
arms of the galaxy. In this section we describe how we model 
both types of field, and show how the random and ordered 
fields are linked by the differential rotation of the galaxy. 

Simulations of a turbulent MHD fluid suggest that 
random magnetic fields can be generated within such fiu- 
ids via turbulence (|Cho et al. 20091 ICho fc Vishniac 20001 
IKulsrud et al. 19 97). These simulations have shown that 
in such a turbulent MHD fluid, equality between random 
magnetic field energy and turbulent fluid kinetic energy is 
reached after several turnover times. This result forms the 
basis for our model of random magnetic fields in galaxies. 
We then model the formation of the large-scale ordered field 
by a simplified model based on the differential rotation of 
the galaxy. 

The mechanism by which turbulence generates a mi- 
crogauss random magnetic field is field line stretching. The 
turbulent fluid motions draw out field lines, and thus, in 
ideal MHD, the field strength along that direction must in- 
crease. This increases magnetic pressure, which opposes the 
stretching process. On a timescale of several eddy turnover 
times, this build up in magnetic pressure suppresses further 
field line stretching, and the process ceases when the random 
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magnetic field energy is equal to the energy in turbulence of 
the fluid (Cho & Vishniac 2000). 

Whilst the random field can be accounted for through 
turbulence in the ISM, a mechanism for generating the or- 
dered field must invoke the differential rotation of the galaxy. 
The galactic dynamo theory is commonly used to explain 
the presence of large-scale coherent ordered magnetic fields 
in spiral galaxies (e.g. Ruzmaikin et al. 1988). In dynamo 
theory, a radial magnetic field is wound into a toroidal field, 
and amplified, via the differential rotation of the galaxy. In 
addition to this so-called "w- effect " , a further process ex- 
ists whereby a toroidal field is transformed into a radial 
field - this is known as the "a-effect" ((Parker 1955). Ar- 
shakian et al. (2009) have presented detailed calculations of 
these processes. In the present paper, a simpler model for 
the generation of an ordered field is used. The generation 
of magnetic fields is broken up into a two step process. In 
the first step, the turbulence in the ISM stretches and am- 
plifies the random magnetic field. In the second step, differ- 
ential rotation winds the random field to create an ordered 
field. Such a division is physically motivated by the relevant 
timescales for the two steps. The random field is amplified 
on a timescale comparable to the eddie turnover time. This 
scale is much shorter than the characteristic timescale for 
ordered field generation, the galactic rotational period. In 
other words, it is a good approximation to treat the ran- 
dom field energy as amplified to a point at which it is equal 
to half the energy in turbulence, before allowing this ran- 
dom field to be converted into an ordered field by the large 
scale differential rotation. Such an approach has the advan- 
tage of being simple while capturing the relevant physics, 
and therefore can be implemented in a self-consistent galaxy 
formation and evolution model. 

The ordering of the magnetic field, performed by the 
mean-scale dynamo, is a much slower process. Arshakian 
et al. find the relevant timescale to be proportional to disk 
size. By contrast, in our approximation the e-folding time 
for this process is simply related to the angular rotation ve- 
locity of the galaxy, and is thus constant at a given redshift. 
It is worth pointing out that the two timescales are simi- 
lar: Arshakian et al. (their Equation 7) give the ordering 
timescale t* = ^ = (-) where the ratio between disk 
scale height and turbulence length scale is - ~ 5, and an- 
gular velocity fi of the galaxy can be expressed in terms of 
the rotation period trot. Therefore, their t* is very close to 
our assumed rotation timescale t TO t- 



2.2.1 Random fields 

As discussed above, since the turnover timescale of turbu- 
lent eddies is much less than cosmological timescales, the 
equality between random magnetic field energy and turbu- 
lent kinetic energy is assumed to hold instantaneously, at the 
point of energy injection. In the absence of an ordered field, 
the instantaneous equality between random field energy and 
turbulent energy demands half the energy injected into the 
ISM goes into turbulent fluid motion, and the other half into 
random magnetic field energy. The processes that inject and 
remove energy from the ISM will now be considered. 

One of the most important sources of energy injection 
into the ISM are supernovae. They inject energy at a rate 



^SFesN, where ^sf is the star formation rate in Moyr -1 
and esN is the energy released by supernovae per solar mass 
of star formation. The quantity esN can be estimated as fol- 
lows. Assuming the minimum stellar mass to be O.IMq, and 
taking all stars with masses greater than 8Mq to terminate 
as supernovae, the supernova rate Rsn can be related to the 
star formation rate by 



-Rsn = *sf (t) 



j^l (p(M)MdM 



*sf(*) Js 



3 dM 



J o °^ M~ 2 - 3B MdM ' 



(1) 



where we have assumed a Salpeter initial mass function 
(j>{M) oc M -2,35 . 

If each supernova injects energy Ek into the ISM, 



esN = 



RsnEk 

*SF 



(2) 



If Ek ~ 10 44 J, approximately the kinetic energy of the 
supernova envelope, then esN ~ 10 41 JMg -1 . 

Star formation removes turbulent energy. If -Eturb is the 
turbulent ISM energy, and M co id the mass of cold gas in the 
galaxy, the turbulent energy removed by the formation of a 
star of mass M* is A * f * -E turo (i.e. the fraction of total gas 
mass removed by star formation, multiplied by the turbulent 
energy). So the rate at which energy is removed by star 



-E tu 



formation is 

Gas accreting from the dark matter halo deposits its po- 
tential energy in turbulence. If a gas packet of mass AM coo i 
falls onto the galactic disk from a radius r, it adds energy 
j Tt GM ( <r i) AJVI cooi ^ r ^ the cold gas in the disk; we will as- 
sume all this energy goes into turbulence. If we sum over 
all such gas packets, and then differentiate with respect to 
time, we find that the energy injection rate from accretion 

• r> ( m GM{<r)M, 

is E; (/ £ 



ss^dr). 



The total rate at which energy is injected into turbu- 
lence is equal to the sum of the three rates above. In front of 
each term in the equation we add factors which reflect the 
efficiency with which energy is added or removed from the 
ISM by a particular process. For the accretion this factor 
is labelled e gr av, for supernovae feedback esnc and for star 
formation esf. Half of the energy injected into the ISM goes 
into turbulent motions of the fluid, and the other half into 
the random magnetic field energy E™ , so the equation for 
the evolution of the random field energy in the absence of 
differential galactic rotation is 



dt 



esNc^SFesN 



-Eturb , 
eSFTT *SF 



+ e s 



A/cold 

GM(< r)M coo i,i 



dr 



(3) 



where the sum is over the gas parcels cooling at the timestep 
of interest. Here, $sf, M co id and M coo i are calculated self- 
consistently from the galaxy formation model. Note that in 
the absence of rotation (i.e. when the whole field is in the 
random component) we have — tt- 
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2.2.2 Ordered fields 



3 PRELIMINARY RESULTS 



As discussed above, we assume that the ordered field is cre- 
ated through galactic differential rotation. If at time t = 
the energy in the random field is E§jj, assuming that there 
is no further energy injection, the energy in the ordered field 
at time t, Eg d (t), is 



Eg d (t) = E. 



ran 




~t/T old 



(4) 



where the timescale for energy transfer from the random 
to ordered field is r or d. We expect r or d to be comparable with 
t rot , the rotational period of the galaxy, and therefore we set 

Tord — trot • 

In reality, energy is constantly injected into the ISM 
(by the processes considered above). Therefore we must ex- 
tend the above calculation to determine how the energy in 
the ordered magnetic field varies when energy is continu- 
ally injected to the random magnetic field. Let the ran- 
dom field energy be E^ n (t) at time i. At t + At the en- 
ergy is Eg a (t + At) « Eg n (t) + ( dE f°, (t,) ) At . Thus 

( ~~^F^~ " ) ^ ^ as b een added to the random field energy 
in the infinitesimal time interval At. For ( dEB dt } t - j < 0, 
the ordered field is assumed not to change - if the random 
energy declines there is no physical reason for the ordered 
field energy to decline also. If ( —^r- - ^ > 0, using Equa- 
tionUl the ordered field energy at time to due to the random 
field energy added in the time interval [t, t + At] is, 



AK rd (to) = ( **3LW1 ) At (l - e (*"*o)/W ) . (5) 



V dt' 



The total energy in the ordered field can now be ob- 
tained by adding all contributions, 



in., ! , - / ( \, ) f(tMdt, (6) 



where, 



/ (t/trot) = 



1 - e 



(t-t )/*rot 



for^>0 
dt 

for ^ < 
dt ~ 



, (7) 



Equation [3] for the random field energy holds with an 
extra loss term included to reflect the fact that random field 
energy is continuously being converted into ordered field en- 
ergy. The final equations for the evolution of the magnetic 
field energies are 



dt 



dt 



dt 



1 rf-fiturb 

2 dt 



-/(*/** 



dEZ 



dt 



(8) 



(9) 



There are four free parameters in the model. Three param- 
eters describe the efficiency with which turbulence is added 
to/removed from the cold gas by gravitational infall (e gra v), 
supernovae feedback (esNc) and star formation (esf)- These 
determine the total magnetic field strength. The fourth pa- 
rameter is a timescale relating to the rate at which the ran- 
dom magnetic field is transferred to the ordered field, r or d, 
which in the previous section we have set to the rotation 
period of the galaxy. 

We set esF = 1, as physically one would expect all tur- 
bulent energy associated with the cold gas to be removed 
once the gas turns into stars. Figure [T] shows the interplay 
between the three mechanisms with this value fixed for a 
range of halo masses. Shabala & Alexander (2009) estimate 
that about 1 per cent of SNe energy must go towards heating 
and ejecting galactic gas in order to reproduce the low-mass 
end of the present day stellar mass function. We therefore 
expect esNc to be of this order, i.e. esn c ~ 0.01. Figure [T] 
shows that in this case gravitational infall and star forma- 
tion dominate supernovae feedback, and thus the exact value 
of esNc is unimportant. We thus set esNc = 0.01. 

As halo mass increases, so do the contributions to tur- 
bulence from gas cooling and the removal of turbulence due 
to star formation. In low mass haloes (Mhaio < 10 12 M©) at 
the present epoch, as halo mass increases, more turbulence 
is added through gas cooling than is removed by star forma- 
tion. While significant amounts of gas can cool, an increase 
in the size of the galaxy disk means the gas surface density 
drops, and thus star formation is less efficient. As a result, 
the magnetic field strength increases with mass (Figure [2J. 

The turbulent energy density and magnetic field 
strength peak at Mhaio ~ 10 13 Mq. At these masses, both 
supernovae and AGN feedback are too weak to quench 
large amounts of gas cooling. Massive haloes are more likely 
to host powerful AGNs (Shabala et al. 2008), and thus in 
higher mass haloes (Figure [TJ> and h) AGN feedback trun- 
cates gas cooling and star formation at early epochs, with 
the turbulent energy density remaining unchanged since 
z ~ 1. Because the gas cooling time increases with decreas- 
ing density, powerful AGNs in massive hosts play a more 
important role in quenching gas cooling and star forma- 
tion than do their lower power counterparts in less massive 
haloes. In the most massive haloes (M ha io > 10 15 M ) AGN 
feedback cannot completely quench gas cooling, and thus the 
magnetic field strength begins to increase once again. 

The average magnetic field strength B in a galactic disk 
is a function of the turbulent energy density. Thus, apart 
from the amount of turbulence present, another important 
parameter is the size of the disk. In the following section 
we use observations of local galaxies to constrain disk sizes, 
and consequently find the value of egrav ~ 0.003. This is the 
value adopted for the rest of the paper. 

Figure [3] shows the evolution of the magnetic field 
for each of these haloes. The volume of the radio disk 
used to calculate magnetic field strengths was taken as 
7ri?diakTdisk, where the disk radius is i?di s k = 0.1i? v ir 
([Shabala fc Alex ander 2009 ) and we adopt disk thickness of 
Tdisk = 0.6 kpc (|Fitt fc A lexande r 1993|l . In practice, disk 
radii derived from radio observations are subject to selection 
effects. We return to this point in the following section. Key 
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Figure 1. Total turbulent energy in the cold gas. Left panels: e grav = 0.003. Right panels: e grav = 0.001. All curves are for esNe = 0.01 
and esF = 1- Star formation and cold gas infall dominate the SNe feedback contribution to turbulence. In high mass haloes, early 
truncation (at z ~ 1) of gas cooling, followed after a time lag by a truncation in star formation, results in the spike in, and the eventual 
flattening of, the turbulent energy curve. 
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Figure 2. Magnetic field strength at the present epoch as a function of halo mass for e gra v = 0.003 (solid line), e gra v = 0.01 (dashed) 
and e g rav = 0.001 (dotted line). All models have egp = 1, esNc = 0.01. 
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Figure 3. Contributions to the total magnetic field (solid lines) from the random (dotted) and ordered (dashed) components as a 
function of halo mass. All models have egp = 1, esNo = 0.01, e grav = 0.003. 



features of the model are evident in Figure [3] In low-mass 
haloes, interplay between the infall of cool gas and star for- 
mation result in substantial contributions to the turbulent 
energy budget right up to z — (Figure [T|. As a result, 
the bulk of the turbulent (and hence magnetic field) energy 
density is in the random component. As halo mass increases, 
so does the strength of AGN feedback. Consequently, the 



most massive galaxies (Figure Oi) are red and dead, with 
the bulk of the cooling and star formation taking place at 
z > 1. Thus, by z = a greater fraction of the turbulent 
energy has had time to be transferred to the ordered field 
than in their lower-mass counterparts. 
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4 MAGNETIC FIELDS IN GALAXY DISKS 
4.1 Disk sizes 

So far we have assumed that radio emission is uniform across 
the disk. Contrary to the assumptions used to plot Figures[3] 
and [21 in real disks the emissivity typically decays exponen- 
tially away from the disk centre (Freeman 1970). As a re- 
sult, essentially all the source flux is observed within a few 
scale heights, while the observed radial extent of the disk 
will be greatly affected by selection effects such as the limit- 
ing beam-averaged surface brightness. In general, it is diffi- 
cult to quantify these effects, and we choose to parametrise 
the size of the radio disk via the quantity edisk, such that 

-Rdisk, radio — ^disk-^disk • 

Figure 0] compares predicted disk sizes with observa- 
tions of late-type galaxies from the sample of Fitt & Alexan- 
der (1993). These authors considered radio properties of 
an optically-complete subsample of 165 late-type galaxies 
brighter than magnitude Bt = +12, drawn from the Re- 
vised Shapley Ames Catalogue (RSA; Sandage & Tammann 
1981). These were complemented by VLA observations at 
1.49 GHz (Condon 1987), yielding an optically-complete 
sample of 146 galaxies with resolved radio disks. Of the 
remaining 19 galaxies, one had unresolved radio stucture, 
eight more had no radio detections; and the remainder suf- 
fered from confusion or were not properly imaged by the 
observations. 

Stellar masses in Figure [4] and the subsequent discus- 
sion were derived from A'-band magnitudes under the as- 
sumption that all stars are of solar type (Shabala et al. 20081 
Nikol icTet al. 2004)1 . Comparison of predicted and observed 
disk sizes shows that setting edisk = 0.25 gives good agree- 
ment with observations. It is worth noting that since smaller 
disks are associated with higher turbulent energy densities 
(for fixed total turbulent energy), greater magnetic field 
strengths are predicted for these. 



4.2 Magnetic field strengths 

Fitt & Alexander (1993) used the observed angular sizes 
of the radio disk to derive equipartition magnetic field 
strengths. They modelled the synchrotron disk in each 
galaxy as having the radius given by imaging the galaxy 
at 1.49 GHz; and took an equivalent disk width of 0.6 kpc. 
This is the disk thickness we adopt for the rest of this paper. 
With disk volume fixed, magnetic field strength is a function 
of a single parameter e gra v, parametrising the efficiency with 
which the potential energy of infalling cold gas is converted 
to turbulent energy. The predicted and observed magnetic 
field strengths are shown as a function of stellar mass for 
different values of e gra v in Figure [5] 

VLA sensitivity sets a limit on the strength of the weak- 
est detectable magnetic field in the Fitt & Alexander (1993) 
sample. This is given by (see Fitt & Alexander Equation 2) 



* in" 1 



B n 



MG 



384 



f #fwhm/2 \ 
Jy \ arcsec J 



2/7 



(10) 



where the coefficient is different to the Fitt & Alexander 
value due to their assumption of cosmic ray energy densities 
being negligible (ry p ~ 0). By contrast, we assume a ratio be- 
tween proton and electron energy densities of r/p = 100 that 
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Figure 6. Distribution of total magnetic field strengths. Top 
panel: e gra v = 0.003. Bottom panel: e grav = 0.01. Data points 
are again from Fitt & Alexander (1993). The vertical solid line 
shows the minimum field strength required for the object to have 
detectable radio emission. Observed points below this value corre- 
spond to RSA galaxies with no radio identifications. Model curves 
show the mean predicted value (solid line) and upper and lower 
limits (dashed lines) due to uncertainties in the observed stellar 
mass function. 



is typical for our Galaxy (e.g. Longair 1994 and references 
therein), yielding an extra factor of 3.7. Beam Full Width at 
Half-Maximum is #fwhm ~ 1 arcmin, and limiting surface 
brightness E m i n = 0.35 mjy. This yields -B m in = 5.7 /xG. 

Figure [6] plots the distribution of magnetic field 
strengths. These were obtained by binning the Fitt & 
Alexander galaxies in stellar mass, and convolving this dis- 
tribution with the results of Figured Our model (left panel) 
provides a good match to the data at low and intermediate 
magnetic field strengths, corresponding to low and interme- 
diate mass galaxies. The model overpredicts the counts at 
the bright end, primarily due to the expected peak in B at 
2 x 10 11 M Q (Figure(S]|. A model with e gra v a factor of 
3 higher than in our best fit model is included for compari- 
son (right panel) . Counts at low values of B are significantly 
underpredicted, consistent with this model predicting higher 
values of B at a given stellar mass than found by Fitt & 
Alexander (Figure [5]). 
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Figure 5. Distribution of total magnetic field strengths 8 as a function of stellar mass. Curves are for models with e gra v = 0.001 
(dotted), 0.003 (solid) and 0.01 (dashed). Data points are from Fitt & Alexander (1993) corrected for a non-negligible proton energy 
density (see text). Radio disk size is fixed at i?disk radio = 0.25-Rdisk as given by Figuref4] All models have egp = 1, egNe = 0.01. 



5 LOCAL RADIO LUMINOSITY FUNCTION 
IN LATE- TYPE GALAXIES 

5.1 Synchrotron luminosity 

Synchrotron luminosity is the most straightforward obser- 
vational manifestation of magnetic fields. Following Longair 
(1994), assuming equipartition this luminosity is given by 



^synch 



2 B 7/2 V 
3/xo G{a)ri P 



(11) 



where V is the volume of the emitting region and r] p the 
ratio between proton and electron energy densities. Con- 
sistent with observations of our Galaxy (Longair 1994 and 
references therein) we adopt rj p = 100. The constant G(a) 
depends weakly on spectral index a, and the minimum and 
maximum energy cutoffs v m i n and ^Wx- F° r an electron 
power-law distribution N(E) oc E~ p where p — 2.5, we have 
a = (p — l)/2 = 0.75. Adopting i/ m i n = 10 MHz, we have 
(Longair 1994) 



G(a) = 4.15 x 10 1 



5.2 Radio luminosities 



GHz 



(12) 



The Fitt & Alexander (1993) sample used in the preceding 
section is representative, but by no means complete. Using 
radio luminosity as a proxy for magnetic field strength allows 
us to use a complete (in both optical and radio) sample of 
Shabala et al. (2008). Synchrotron luminosity for late-type 
galaxies in the sample is plotted as a function of stellar mass 
in Figure [7] Here, late-type galaxies are defined as galax- 
ies with concentration indices in the Petrosian r-band of 
C > 0.375. Stellar masses were obtained from z-band mag- 
nitudes by assuming solar-type stars (Shabala et al. 20081 
INikolic et al. 2004[> . 

Given a disk size (i.e. the parameter tdisk, set to 0.25 
in Section 14. ip , the model predicts a single radio luminos- 
ity for each halo, and hence a given stellar mass. Figure [7] 
shows that these agree well with mean observed luminosities 
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Figure 7. Radio luminosity as a function of stellar mass. Top 
panel: All disk-type galaxies with radio detections. Data are from 
Fitt Sz Alexander (1993; filled squares) and Shabala et al. (2008; 
dots). The Fitt & Alexander observations at 1.49 GHz were con- 
verted to 1.4 GHz values assuming a spectral index of 0.8. Bottom 
panel: Mean radio luminosity. Filled squares are again for the Fitt 
& Alexander sample; open circles denote the Shabala et al. (2008) 
sample. The Shabala et al. sample is flux- and volume-limited, 
thus V/V max corrections must be applied. The resultant corrected 
mean luminosities are shown by open squares. The minimum de- 
tectable luminosity at the front of the volume (~ 6 X 10 21 W Hz -1 
at 2 = 0.1) is shown by a horizontal black line. No such cor- 
rections are required for the surface-brightness limited sample 
of Fitt & Alexander (1993). Curves shows model predictions for 
£grav = 0.003 (solid), 0.001 (dotted) and 0.0003 (dashed). 



at low stellar masses. However, at higher masses the radio 
luminosity is overpredicted. The radio-loud galaxies in the 
flux-limited sample of Shabala et al. make up the bright tail 
of the radio luminosity distribution. The majority of galax- 
ies have radio luminosities below the detection limit, and 
thus the mean radio luminosity at a given stellar mass was 
calculated by assuming a Gaussian distribution with the ob- 
served scatter of 0.3 dex. By constrast, the Fitt & Alexander 
observations were performed with much greater sensitivity, 
and thus no such corrections are necessary. 

The peak radio luminosity in Figure [Jj corresponds to 
the B-field peak in Figure [21 and is related to an increase in 
gas cooling (and hence the random component of the mag- 
netic field) at late times in haloes with Mhaio ~ 10 13 Mq 
(Figures QJ and f, and Figure [3J:). This late cooling is a re- 
sult of the way AGN feedback is implemented in our model. 
In haloes more massive that ~ 10 12 Mq supernovae feedback 



is too weak to quench gas cooling and star formation. AGN 
heating can in principle do this, however only the dense cen- 
tral regions can be heated by weak AGNs found in haloes 
with mass < 10 13 M . These cores have short cooling times, 
and thus the net effect of AGN heating is negligible. By con- 
trast, in the most massive hosts powerful AGNs can heat rar- 
efied gas at larger radii, where the cooling times are longer, 
and thus AGN heating is more effective. 

As pointed out by Shabala & Alexander (2009), an im- 
portant feedback ingredient missing in our model is gas ejec- 
tion by powerful radio sources. In addition to shock heating, 
jet-inflated radio cocoons sweep up the ambient gas through 
which they expand, transporting it to larger radii (e.g. 
Alexander 2002; Basson & Alexander 2003). The longer cool- 
ing times will then prevent the hot gas from cooling, decreas- 
ing the turbulent energy densities, magnetic field strengths 
and radio luminosities in AGN hosts. A comparison of the 
present day cold gas mass function predicted by our model 
with the observed H I mass function (Z waan et al. 2003[) in- 
dicates that the amount of cold gas in massive hosts is over- 
predicted by 0.5 — 1 dex. The effects of adding AGN gas 
uplifting to our model on the predicted magnetic fields can 
be estimated by decreasing the parameter e gr av by this fac- 
tor. Figure [Jj shows that this provides a better fit to ob- 
servations, however the rise in radio luminosity between 
M* = 10 12 - 3 x 10 12 M is still too steep. This is not sur- 
prising, since this is precisely the range where distributed 
supernovae feedback cannot stop large amounts of gas from 
cooling, while AGN heating is not yet powerful enough to 
affect gas at large radii with long cooling times. One would 
expect the gas uplifting to play a more important role in 
these intermediate-mass haloes. 



5.3 Radio luminosity function 

The stellar mass function for local late-type galaxies was 
constructed from the Shabala et al. sample. This stellar mass 
function was then converted into a halo mass function for 
late-type galaxies by using the best-fit model of Shabala & 
Alexander (2009). The radio luminosity function can then 
be constructed from the halo mass function and predictions 
for individual haloes (Figure [Jj. 

The resultant RLF is shown in Figure [8] As in Figure [Jj 
the number of bright sources is vastly overpredicted. De- 
creasing the amount of cold gas by a factor ~ 10 reconciles 
the model with observations, suggesting outward gas trans- 
port by powerful radio sources is important to the evolution 
of galactic magnetic fields. Observationally it should be pos- 
sible to distinguish between the gas heating versus ejection 
scenarios through use of X-ray data. Frequent AGN heating 
and subsequent cooling (due to short cooling times) in dense 
central regions would manifest itself in copious amounts of 
X-ray emission. By contrast, gas transport to large radii 
would simply suppress cooling, and thus not give rise to any 
appreciable X-ray signatures. 



6 SUMMARY 

We have developed an analytical model to follow the cosmo- 
logical evolution of magnetic fields in disk galaxies. This is 
done by assuming equality between the magnetic field energy 
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Figure 8. Radio luminosity function for late-type galaxies. Model 
predictions are again given for e gra v = 0.003 (solid curve), 0.001 
(dotted) and 0.0003 (dashed). Observed points are from Fitt & 
Alexander (1993; filled squares) and Shabala et al. (2008; open 
circles). V/V max corrections are once again applied to the Shabala 
et al. sample. 

density, and the turbulent energy density in cold disk gas. 
The gas cooling and star formation histories are followed us- 
ing the Shabala & Alexander (2009) galaxy evolution model, 
which includes a physically-motivated prescription for AGN 
feedback, together with supernovae and reionization feed- 
back. The model successfully reproduces the observed stellar 
mass functions and shutdown of star formation in massive 
galaxies out to redshifts of 1.5. 

Three mechanisms alter the turbulent energy budget in 
the cold gas. Star formation removes turbulence associated 
with gas parcels that collapse to form the stars. The turbu- 
lence is increased by gravitational infall of cool gas onto the 
disk; and also by high-mass stars driving shocks through 
the ISM as they end their lives in supernovae explosions. 
Two types of magnetic fields are identified: the random field, 
which is incoherent on scales comparable with the turbulent 
eddy scale; and the ordered field, generated from the random 
field by differential rotation of the galaxy. 

Two local samples are used to test the models. The 
model reproduces magnetic field strengths and radio lumi- 
nosities well across a wide range of low and intermediate- 
mass galaxies. However, the radio luminosities and magnetic 
field strenths are overpredicted significantly in high mass 
galaxies due to an overprediction in the amount of gas cool- 
ing. Inclusion of outward gas transport by powerful radio 
sources is required to reconcile the model with observations. 
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